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A computational region of arbitrary cross section presents a significant 
problem in the generation of a mesh. Simple orthogonal meshes are 
difficult to use because the mesh points do not naturally fall on the 
region’s boundaries. Differencing and interpolation schemes become 
complex and cumbersome, and it is difficult to extend these schemes to 
higher order because of the complex logic required. Higher order schemes 
are desirable as they allow calculation of a flow to a given level of 
accuracy with a lower mesh density and hence less storage than a lower 
order scheme. High accuracy solutions are possible for a region of 
arbitrary cross section when a boundary-fitted computational mesh is 
employed. A boundary-fitted mesh is defined as a mesh in which the 
boundary (i.e., a duct wall) is coincident with the mesh points that are 
used for finite difference expressions at, and adjacent to, the 
boundary. Interpolation is not required, and extension to higher order 
differencing is straightforward. This is a significant benefit when the 
boundary conditions have a dominant influence on the solution. 


This paper will discuss the application of Smith and Wiegel's method for 
generating boundary fitted coordinate systems (discussed in their 
AIAA-80-0192 paper entitled, "Analytic and Approximate Boundary Fitted 
Coordinate Systems for Fluid Flow Simulation") for two practical flow 
problems characterized by complex surface geometry: 
o radial mixer lobe 

o subsonic inlet designed for high angle-of-attack capability 
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Figure 4.- 


Inlet contour. 


335 


in the method of Smith and Wiegel, two disconnected boundaries are 
defined and an explicit functional relation is used to establish the 
transformation between the physical domain and the computational domain. 


The physical domain is defined by a cartesian coordinate system; the 
computational domain is defined with the variables £ , n and C with the 
values: . _ 

^ l 

Oi T] t i. 

C ‘ 1 

Two possible connecting functions are suggested: linear and a cubic 

parametric polynomial. The following cubic polynomial equation was used 
to generate meshes for both the lobe mixer and the subsonic inlet: 
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where : 


VVi )>Y / <Vi >- Z VV^ >• x ■ >’ 2 » re the 

boundary points in the physical domain 
dX< (V^). «/ (VI > 4 *< <VU ./■ 1.2 are the 

d^ df\ 

derivatives of the boundary points in the physical domain 


fj(n) » 2n 3 - 1 

f 2 (n) - -2n 3 « In' 
f 3 tn) • n 3 - 2n 2 + n 
f 4 tn) - n 3 - n 2 


The cubic connecting function forces orthogonality at the boundaries of 
the physical domain by calculating the derivatives <w 

dY/ (^ .^ ) and dZ/ ^ ^ ^ from the cross product of the tangential 
derivatives and then dividing by the magnitude of the normal vector. 


Four extensions of the Smith and Wiegel method were necessary in order to 
successfully apply their technique to the mixer lobe and subsonic inlet. 


First, because of the nature of the mixer and inlet geometries, points 
defining the boundaries had to be positioned using a geometric 
progression. 


S = a + ar + ar 


. + ar 


N-l 


«(!-/) 

1-r 


where 

S = the total length of the boundary 
a = first increment 
r = scale factor 

N = number of cells (one less the number of boundary points) 




For the mixer, the scale factor r was varied linearly from r - 1 at the 
mixer entrance plane (where the boundary is an arc) to r = r max a t 
mixer exit plane (where the boundary is highly distorted). This makes it 
possible to force the mesh to migrate to regions of interest without 
causing significant distortions in the mesh from plane to plane. The 
optimal distribution of mesh occurred when the upper and lower boundary 
mesh points were stretched in opposite directions. 
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SECONDARY LOBE - PLANE • 94 

GEOMETRIC PROGRESSION LOWER BOUNDARY 1.0000 UPPER BOUNDARY 1.0000 
SLOPE SCALER LOWER BOUNDARY 1.0000 UPPER BOUNDARY 1.0000 
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SECONDARY LOBE - PLANE J 07 

GEOMETRIC PROGRESSION LOWER BOUNDARY 0.9750 UPPER BOUNDARY 1.0500 
SLOPE SCALER LOWER BOUNDARY 1.0250 UPPER BOUNDARY 1.0500 



SECONDARY LOBE - PLANE 120 

GEOMETRIC PROGRESSION LOWER BOUNDARY 0.9500 UPPER BOUNDARY 1.1000 
SLOPE SCALER LOWER BOUNDARY 1.0500 UPPER BOUNDARY 1.1000 

Figure 5.- Geometric progression for boundary points for secondary stream. 
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PRIMARY LOBE - PLANE 94 

GEOMETRIC PROGRESSION LOWER BOUNORRY 1.0000 UPPER BOUNORRY 1.0000 
SLOPE SCALER LOWER BOUNDARY 1.0000 UPPER BOUNDARY 1.0000 



PRIMARY LOBE - PLRNE 107 

GEOMETRIC PROGRESSION LOWER BOUNORRY 1.0500 UPPER BOUNORRY 1.0500 
SLOPE SCRLER LOWER BOUNORRY 1.0500 UPPER BOUNORRY 1.0500 


ORIGINAL PA JR IS 
OF POOR QUALITY- 



PRIMARY LOBE - PLRNE 120 

GEOMETRIC PROGRESSION LOWER BOUNORRY 1.1000 UPPER BOUNORRY 1.1000 
SLOPE SCALER LOWER BOUNORRY 1.1000 UPPER BOUNDARY 1.1000 

Figure 6.- Geometric progression of boundary points for primary stream. 
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The inlet has certain regions (hilite, throat, etc.) which require a fine 
computational mesh to insure a detailed analysis. For this reason, four 
regions along each inlet contour and five regions along the boundary of 
the analysis domain required individual geometric progressions. The 
scale factor, r , and the number of cells, N, of each region must be 
chosen to insure a smooth progression in cell length along each of the 
boundaries. 



Figure 7.- Geometric progression regions along 
inlet contour and analysis boundary. 
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The second extension uses a ramping function to regulate the dependence 
of the connecting function on the boundary slope. This connecting 
function is an explicit functional relation used to establish the 
transformation between the physical domain and the computational domain. 

For the mixer lobe, this dependence was regulated to redistribute the 
internal mesh points and reduce mesh skewness. 

In the case of the subsonic inlet, it was found that a constant value for 
each plane was sufficient to insure against mesh line cross-over. 



Without ramping function 



With ramping function 

Figure 8.- Connecting function dependency 
on boundary slope. 
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The third extension utilizes the concentration function suggested by 
Smith and Wiegel, but uses it to force the mesh in the direction of both 
boundaries of the mixer lobe. More mesh was then needed to be linearly 
added to fill the void created by this mesh concentration. 

The inlet only required the mesh to be forced towards the inlet contour. 

A concentrated mesh was assumed unnecessary along the spinner boundary; 
it was felt that for a potential flow analysis the flow about the spinner 
would not propagate upstream and affect the solution at the regions of 
interest (hilite, throat, etc.). The mesh concentration for both the 
mixer and the inlet permits flow analysis within the boundary layers. 



Mesh concentrated towards 
inner boundary 



Mesh concentrated towards 
outer boundary 



Mesh concentrated towards 
both boundaries 

Figure 10.- Mesh concentration. 



The fourth extension applies to the subsonic inlet only. It was 
necessary to produce a computational mesh which possessed a smooth 
progression of cell metrics and cell volumes in all directions to allow a 
solution process of a flow analyser to use the grid efficiently. The 
interior points of the computational mesh were "smoothed" by a multiple 
application of a five point diffusion operator: 


X(L,l) new = CL |x(L * 1,1) + X(L + 1,1) +X(L, I- 1) +X(L, I + 1 ) - 4 * X(L, I) 0 id) 

Y(L, l) new = a{Y(L-1,l) + Y{L+1,l)+Y(L,l + 1) + Y (L, I + 1} - 4 . Y ( L. I) old ] 


The value of 01 and the number of times of application were determined by 
trial and error. 

The "smoothed" boundary points could not be determined from the five 
point diffusion operator since one of the required smoothing points would 
be outside the mesh region. Their values were determined from the 
intersection of, the lines defined by the "smoothed" interior mesh points 
and the boundaries. 


X{U) new - Cl{x(L-l,,) + X(L 

YIL - l} new - a [y(L- 1, /} + y (L 


+ U) +X(L, I- i) +X( L( 
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Figure 12.- Five 


Poi'nt di ff us1on operator. 



a = 0.05; Number of interations = 1 



a = 0.05; Number of interations = 3 
Figure 13.- "Smoothed" computational mesh. 







a = 0.05; Number of interations = 12 
Figure 13.- Concluded. 
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slope of segment (?) 

Y(L, I) - Y(LJ-fl) 

I) - X(L f I + 1) 

equation of segment © 

Y-Y(L,n = Mifx-X(L,l)J © 

slope segment © 

Y(L + 1, I + 1) - Y(L + 2, 1 +1) 

M 2 = X(L+1,I + 1)-X(L + 2,I + 1) 

equation of segment © 

Y-Y(L + 1,l+1) = M 2 £x-X(L + 1,l + 1)J © 

since a line thru segment © intersects segment © r 
the X's and Y's of equations © & © equal each other. 

Solving for X: . * 

X(L. I ♦ 1)naw° Ml [ X,L - "j* 2 1 XI ' 1 * * T Y(L + ^ 1 + 11 • Y(l - 11 

Ml - Mi 

Solving for Y: 

Y(L, I + 1) new = M 1 [ X(L ' 1 + ^new' X < L ' »] + Y(L, I) 
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SECONDARY LOBE - PLANE J20 

GEOMETRIC PROGRESSION LOWER BOUNDARY 0.9500 UPPER BOUNDRRY 1.1000 
SLOPE SCALER LOWER BOUNDARY 1.0500 UPPER BOUNDARY 1.1000 




PRIMARY LOBE - PLANE 120 

GEOMETRIC PROGRESSION LOWER BOUNDARY 1.1000 UPPER BOUNOARY 1.1000 
SLOPE SCALER LOWER BOUNDARY 1.1000 UPPER BOUNDARY 1.1000 


Figure 14.- Example mesh for last mixer plane. 
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Figure 15.- Example mesh for 
subsonic inlet. 


Conclusions 


The method of Smith and Wiegel can be used to generate meshes for mixer 
lobes and subsonic inlets that are compatible with flow analysis codes 
requiring a boundary fitted coordinate system. Successful application of 
this mesh generator required development of procedures to distribute the 
mesh points along the boundaries, to regulate the dependence of the 
connecting function to the local boundary slope, to concentrate the mesh 
into regions of special interest, and to modify the mesh grid so that it 
possessed a smooth progression of cell metrics and cell volumes in all 
directions. The method of Smith and Wiegel when coupled with the 
extensions mentioned above has proven to be easy to use and control for 
the inlet and mixer lobe geometries investigated. 

The next step is the formulation of a truncation error monitor for 
arbitrary meshes. This monitor will define where in an analysis domain 
the grid length scales must be changed and by what amount in order to 
equalize truncation errors over the entire analysis domain. Once these 
errors have been equalized, this same monitor will use several levels of 
grid distribution (of the above analysis grid) to then make estimates of 
the absolute truncation error spectrum. This work is currently under 
contract with the NASA Langley Research Center. 
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